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Abstract. A seminal technique of theoretical physics called Wick's theorem interprets 
the Gaussian matrix integral of the products of the trace of powers of Hermitian matri- 
ces as the number of labelled maps with a given degree sequence, sorted by their Euler 
characteristics. This leads to the map enumeration results analogous to those obtained 
by combinatorial methods. In this paper we show that the enumeration of the graphs 
embeddable on a given 2-dimensional surface (a main research topic of contemporary enu- 
merative combinatorics) can also be formulated as the Gaussian matrix integral of an 
ice-type partition function. Some of the most puzzling conjectures of discrete mathe- 
matics are related to the notion of the cycle double cover. We express the number of 
the graphs with a fixed directed cycle double cover as the Gaussian matrix integral of an 
Ihara-Selberg-type function. 

1. Introduction 

Let M = (Mij) be an N x N Hermitian matrix, i.e., Mjj = Mji for every 1 < i,j ' < N, 
where Mji denotes the complex conjugate of M^. Let <p(M) = J2i a iYlij£i ^ij be a 
polynomial in its entries, where I ranges over a finite system of multisets of elements of 
N x N and a/s are complex numbers. We start with the following notation. 

Notation 1.1. Denote (■) by 

(<p) = J2 ai <U M ^ = E a 'E II ( M * M *>' 

I ijel I P (p,q)&P 

where P ranges over all partitions of I into pairs, and forp = (pi,P2), q = Q2) we. have 
(MpM q ) is non-zero if and only if p± = qi and p 2 = q\ and in that case (M p M q ) = 1/N. 



1.1. The motivation. Wick's theorem interprets (tp) as the Gaussian matrix integral (see 
Section [2]). This has been successfully used in enumeration of combinatorial structures and 
in other fields of mathematics [21 El HI El El El HE1 [IH1 EQl EHl EH [25l EH EI]. 

The basic combinatorial structure studied in this paper is that of a graph. A graph G is 
a pair (V, E) where V is the set of vertices and E is the set of pairs of vertices called edges. 
The vertices of an edge are called its endvertices. A graph may have multiple edges, i.e., 
different edges with the same endvertices, and also loops, i.e., edges whose two endvertices 
are equal. The graphs without multiple edges and loops are called simple. A map is a 
graph together with a fixed cyclic ordering of the incident edges of each vertex: it defines 
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an embedding of the graph on an orientable 2-dimensional surface (see [21]). Note that a 
map is a graph that is embedded on a given 2-dimensional surface. 

Principally studied functions ip are products of powers of the trace of M % for a non- 
negative integer i. For such functions ip, (<p) has a useful graphic interpretation as the 
number of labelled maps with given degree sequence, sorted by their Euler characteristics 
(these maps are the Feynman diagrams for the matrix integral) . This seminal formula was 
first described in [5]. While the Gaussian matrix integral of a polynomial in the traces of 
the powers of M is well understood, the Feynman diagrams expansions of the Gaussian 
matrix integral of the exponential of such polynomials hold in general only on the level 
of formal power series. We will denote this formal equality by ' = '. The expression 
f(x) ' = ' g(x) simply means that all the derivatives of /, g are equal when the variable 
x is set to zero. Similarly the expression f(x) ' < ' g(x) means that all the derivatives 
of f,g satistfy the inequality when the variable x is set to zero. Finally, let s(N) be a 
positive integer function in N and let fj^(x), gi^(x) be two formal power series in variable 
x, whose coefficients depend on a positive integer N. The expression fpf(x) ' — S (N) ' 9n(x) 
(resp. /jv(a^) ' < s (jv) ' 9n{x)) means that all the k-th derivatives of /, g, for k < s(N), 
satistfy the equality (resp. inequality) when the variable x is set to zero. These notations 
are defined analogously for the multivariate generating functions. We describe the classical 
application of the matrix integrals in the next subsection. 

1.2. Enumeration of maps. Let "Hn denote the set of all N x N Hermitian matrices. 
We consider a function ip which maps M G Hn to 

ip(M) :=e <-^ 

where N is the set of all positive integers and z^'s are formal variables. 

Taking the logarithm of (i/j) we get the following formula for connected maps sorted by 
genus as discussed by Brezin, Itzykson, Parisi and Zuber in [5]. 

Theorem 1 (0). 

log w '=' E n2 ~ 29 E n { -^r M ^ ---^ 

9>0 (m,- ,n k )£N k i=l % ' 

where M g (ni, ■ ■ ■ ,nk) denotes the number of connected maps with genus g and vertices 
of degree i for 1 < i < k. 

The convergence of the generating function in Theorem [T] was studied extensively [HI U\ 
EH HH HSJ [22] • For example, the following theorem was proved by Ercolani and McLaughlin 
in [6]. (For more general results, see a lecture note by Guionnet [15] and references therein.) 

Theorem 2 ([6]). Let a function ip map each Hermitian matrix M e 7i N to 
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where I is a finite set of positive integers such that its maximum element is even. Then for 
each i E I there is q > so that for Zi a non-zero real variable with z\ G (0, e*), Theorem 
U\ holds as equalities between analytic functions (in variables Zi) when i/j is replaced by ip. 

1.3. Main contribution. A natural question is whether the Gaussian matrix integral 
method may as well be applied to the enumeration of graphs that are embeddable on a 
given 2-dimensional surface. This enumeration has been a long-standing open problem 
in the enumerative combinatorics, solved only for planar graphs and fairly recently by 
Gimenez and Noy [12] . 

One of our main results (Theorem [3]) is that this enumeration task can as well be formu- 
lated as a Gaussian matrix integral of an ice-type partition function (but whose asymptotic 
analysis is not carried out so far). 

Another main result (Theorem H]) is that the number of graphs with a fixed directed cycle 
double cover is expressed as a Gaussian matrix integral of an Ihara-Selberg-type function. 
Even though the study of the cycle double cover conjectures is extensive, no attempt has 
been made before to calculate these numbers. 

The main contribution of this paper is in applying the matrix integral to combinatorially 
defined functions, in order to loosen the strong connection with maps implied by the 
integration of the functions of traces. In the rest of this section we will show how to 
achieve this goal, introduce necessary concepts, and state our main results more precisely. 

1.4. Directed graphs. A starting idea is to associate directed graphs with matrices. Let 
M G Hn and let D = D(M) = (N, N x N) be a complete directed graph with weights on 
directed edges given by M. 

Definition 1.2. A closed walk in D is defined as a circular sequence of directed edges of D 
so that the tail of each edge equals the head of the consecutive edge. A closed walk is called 
pointed if it has a prescribed beginning. A closed trail is a closed walk which contains each 
directed edge at most once. 

Now Tr(M n ) can be interpreted as Y[ e e P M e , where the sum is over all pointed closed 
walks in D of length n. Similarly, Tr(M 3 ) 4 Tr(M 2 ) 3 can be interpreted as a function which 
maps each Hermitian matrix M G T~Cn to 



where the first sum is over all pointed closed walks p\ in D of length 3, and the second 
sum is over all pointed closed walks p2 in D of length 2. 

Definition 1.3. A subset A of directed edges of D is called eulerian if for each vertex its 
indegree in A is equal to its outdegree in A. We denote by £ = £{D) the collection of all 
eulerian sets in D. 




Pi e€pi 
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Definition 1.4. A subset A of directed edges of D is called even if A can be written as a 
union of edge-disjoint directed cycles of length bigger than two. 

For a set A of directed edges of D, a subset C C A is called a component of A if C 
is a maximal subset of A with respect to inclusion that induces a connected underlying 
undirected graph. 



1.5. Enumeration of embeddable graphs. We say that an integer function f(N) is of 
type 0(N~ l ) if f(N) is bounded above by cN~ l for a constant c and all large enough N. 
Let 7] be a function which maps each Hermain matrix M G 7"^ to 

, w ^,,)=E(n«.) n (Af2) T:r' c ' /2 (1) 

Ae£ eeA C: component of A " 

The first result of the paper is the following expression for planar graphs. 

Theorem 3. Let s(N) < f be any function in N which approaches to oo slower than 
N when N — > oo. The Gaussian integral of the function r/ defined in (QJ) satisfies the 
following. 



z r y n ~'' 
n\ 

n>l r>0 



exp I ^2^2[p{n,r) + d(iV,n,r)] 
e- N \r,(M,Nzy-\y)) 
' < n ' exp ^2 r ) + e2 ( N > n ' r ^ ] 



n! 

,n>l r>0 



where p(n, r) denotes the number of labelled connected simple graphs on n vertices which 
have planar embeddings with r faces, and e\(N,n,r),e2(N,n,r) are functions of type 
O^N- 1 ). 

Let p(n) denote the number of labelled planar connected simple graphs on n vertices. 
Clearly, we have that p{n) = J2 r >oP( n ^ r ) an( ^ t na t p{n)/n\ is equal to the coefficient of 
log<77(M, Nzy~\y))\ z=l . 
It would be very interesting to see if we can recover 

p{n) ~ c n _7//2 7 n n\ 

with c = 0.478 and 7 = 27.2, which is computed by Gimemez and Noy [T2] . 

The embeddable graphs of higher genus contribute to the subleading terms 
ei(N,n,r),e2{N,n,r) in Theorem [31 it should be possible to derive their enumeration 
from our formulas in Section |H 

The function 77 in (P) is an ice-type partition function. Ice-type partition functions play 
a seminal role in modern mathematics and theoretical physics (see [SO]), arid we hope that 
the asymptotic analysis of the matrix integral of the function rj is not out of reach. 
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1.6. Enumeration of directed cycle double covers. Observe that an eulerian set does 
not need to be even. Even sets are closely related to the cycle double cover conjectures. 

Definition 1.5. Let G be an undirected graph. A collection of its cycles is called a cycle 
double cover (CDC) if each edge belongs to exactly two of the cycles. Moreover it is called 
a directed cycle double cover (DCDC) if it is possible to orient the cycles so that they go 
oppositely on each edge. 

Some of the most puzzling conjectures of discrete mathematics are centred around this 
notion. A graph is bridgeless if it cannot be disconnected by deletion of a single edge. 
Clearly a graph with a bridge does not have a CDC. On the other hand, there are 

• Cycle double cover conjecture: Is it true that each bridgeless graph has a CDC? 

• Directed cycle double cover conjecture: Is it true that each bridgeless graph has a 



Our next main result concerns directed cycle double covers. 

Theorem 4. Let A be the set of all pairs (q, K) where q is an even set of edges of D and 
K is a decomposition of q into directed cycles of length at least three. Then 



where the sum is over all isomorphism classes of pairs (G, C) where G is a simple graph 
with at most N vertices and no isolated vertex, and C is a specified DCDC of G. 

We show in Section [5711 that the integrand £ in Theorem H] which maps M e Ti^ to 



is an Ihara-Selberg-type function (see [27], [28]). Thus again, we hope that an analysis of 
the matrix integral of the function £ may be achievable. 

The rest part of the paper is organized as follows. In Section [2] we review the Gaussian 
matrix integral method and Wick's theorem. In Section [3] we introduce identities essential 
to the enumeration of embeddable graphs on a 2-dimensional surface. In Sections H] and [5] 
we prove our main Theorems [3J and HI 



In this section, based on [9], we recall the definition of Gaussian matrix integral, useful 
identities, and Wick's theorem in particular. 
Let M G Hn and let 



DCDC? 





(2) 



2. Gaussian matrix integral and Wick's theorem 



dM = Yl dM u Yl d Re(Afy)d Im(My) 



i i<j 
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denote the standard Haar measure on Hn, where Re (My) and Im(My) denote the real 
part and the imaginary part of M^-. 

For an arbitrary function tp : Hn — > C, the Gaussian matrix integral of p is defined as 

1 



<</?> 



Z (N) 



<p(M)dM, 



(3) 



Hn 



where Z (N) is the normalization factor making (1) = 1, that is, 



Z (N) 



~ NT ^dM. 



Hn 



We are particularly interested in a function p which maps M G Hn to 

</?(M) = J> J] M*;, (4) 

i ijei 

where I ranges over a finite system of multisets of elements of N x N. To compute (</?) 
we introduce the source integral (e 1 *^ 1 ^) for a given N x N Hermitian matrix S, where 
Tr (MS) denotes the trace of the matrix MS. It can easily be computed by 



3 Tr(M5) 



Z (N) 
1 

UN) 

Tr(5 2 ) 

e 2jv 



e -iVIV(^) e lV(MS) dM 



Hi 



-NTr(l(M-§f) e ^P dM 



(5) 

since the trace is linear and Tr(MS) = Tr(SM). On the other hand, for any 1 < i, j < N, 







dS 3l 



D Tv(MS) 



s=o 







dSj, 
d 



Tt(MS)) e^ M ^ 



dS 



j i 



m,n / 



5=0 
Tr(MS) 



5=0 



M, 



i.r 



Thus the derivatives of the source integral becomes 

d d 



dSji dSik 
1 



(e 



Tr (MS) 



Z (N) 
1 



e -iVTr(^-) 



5=0 

d d 



H r 



dS j{ dS, 



■ ■ e 



Tr(MS) 



Ik- 



dM 



5=0 



e- N ^^M l3 M kr --dM 



Zo(N) JHt 
= (MijMkr--), 

where the first equality is due to the Leibniz integral rule. 



(6) 



Using 
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and (15]), we obtain 



(M i:j M kl ■■■) 



© d d 



dSa dSik 



(e 



Tv(MS)\ 



j 3 i 

© d d 



5=0 



dSji dSik 



5=0 



(7) 



and in particular 



{MijMki) 



d d Tr(S 2 ) 

e 2N 



dSji dS^ 







s=o 

2 



d Tr(S 2 )\ Tris^x 



dS^ \dS lk 2N 



e 2JV 



5=0 







9 Y2m,n SmnSi 



dSji \dSik 

e 2iv 



2N 



e 2JV 



5=0 



dSji N 



5=0 



N ' v ' 

Further, it is clear that the derivatives in © and (jHJ) must be taken in pairs (e.g., Sji and 
Sik with I = i and k = j) to get a non-zero contribution. This yields the following result 
known as Wick's theorem. 

Theorem 5 (Wick's theorem). Let M £ T-Ln and I be a multiset of elements of N x N. 
Then 

(U M ^ : 

ijei 



j2 n w«> 

pairing Pel 2 (ij,kl)&P 



n 



N 



(9) 



pairing Pel 2 (ij,kl)€P 

Due to the linearity of the integral and for the function (p defined in 



we have 



l ijei i p (p,q)eP 

where P ranges over all partitions of / into pairs, and for p = (pi,P2), Q = (Qi, Q2) we have 
(MpM q ) is non-zero if and only if p\ = q 2 and p 2 = qi and in that case (M p M q ) = 1/N. 
This is how we arrive at the notation of (•) in Notation ll.il 



3. Trail double cover and pseudoembedding 

In this section we derive several identities, which are essential for applying the Gaussian 
matrix integral method to the enumeration of embeddable graphs on a 2-dimensional sur- 
face. We recall that a map is a graph together with a fixed cyclic ordering of the incident 
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edges of each vertex. This ordering (of each vertex) defines an embedding of the graph on 
an orientable 2-dimensional surface (see [21] )• A map is also called a fat graph; we prefere 
this term since it corresponds to a helpful graphic representation of Wick's theorem. In 
a fat graph F the vertices are made into discs (islands) and connected by fattened edges 
(bridges) prescribed by the cyclic orders. This defines a two-dimensional orientable surface 
with boundary which we also denote by F. Each component of the boundary of F will be 
called a face of F. Each face is an embedded circle (see e.g., [21]). We will denote by G(F) 
the underlying graph of F. 

We denote by e(F),v(F), f(F), c(F), and g(F) the number of the edges, vertices, faces, 
connected components, and genus of F. We recall that 2g(F) = 2c(F)+e(F)—v(F) — f(F). 

In the next sections we will count fat graphs and their relatives. To avoid confusion we 
assume that a fat graph has labelled vertices, i.e., two fat graphs are equal if they are equal 
as sets. We speak about unlabelled fat graphs if the equality is up to isomorphism. 

Let M G H N and we recall that D = D(M) = (N, N x N) is a directed graph with 
weights on directed edges given by M. A closed walk in D is defined as a circular sequence 
of directed edges of D so that the tail of each edge equals the head of the consecutive edge. 
As we have seen in Section IQ1 Tr(M 3 ) 4 Tr(M 2 ) 3 can be interpreted as 

(E n a « 4 (£ n M -r- 

Pi eGpi Vi e&p 2 

where the first sum is over all pointed closed walks p\ in D of length 3, and the second 
sum is over all pointed closed walks P2 in D of length 2. 

A proper pairing of a subset A of directed edges of D is a partition of A into pairs 
ij,ji G N x N of oppositely directed edges. Hence if <p(M) = Tr(M 3 ) 4 Tr(M 2 ) 3 we get 

v> = £ £ II Viv, 

9=91 9a— 97 P (e,e')eP 

where the second sum is over all proper pairings P of the directed edges of the disjoint 
union q = q\q2 ... ^7 of 7 pointed closed walks, from which four have length 3 and remaining 
three have length 2. Two directed edges form a proper pairing if one is reversed the other. 
We also say that such a pair (q, P) contributes to (<p). 

Next we recall that a closed trail is a closed walk which does not repeat edges and 
we denote by VV(r) the set of all subsets of edges of D that can be decomposed into r 
edge-disjoint closed trails, and we let W = U r >o W(r). 

We define a function uo r which maps each Hermitian matrix M G Hn to 

u r {M) := •'/" 2 U XI ( 10 ) 

qeW(r) eG<? 

for each integer r > 1 and L0q{M) = 1. We call y\ g ^ 2 Y[ e eq ^ e the weight of q. 
Next two definitions and proposition are crucial. 

Definition 3.1. Let G be a graph. A trail double cover (TDC) is a collection of closed 
trails so that each edge of G is traversed exactly twice, and in the opposite directions. 
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We remark that each graph has a TDC. Next we introduce a useful notation. A simple 
graph will be called nimple if it has no vertex of degree 0. 

Definition 3.2. Let G be a finite nimple graph with at most N vertices. Then let c(G) 
be the set of all pairs (q, P) so that there is a colouring d of the vertices of G by colours 
{!,..., N}, where each vertex gets a different colour, 



Proposition 3.3. Let r > 0. A pair (q,P) with P a proper pairing of q contributes to 
(u r ) if and only if there is a nimple graph G with a TDC consisting of r closed trails such 
that (q,P) G c{G). 

Proof. If (q, P) G c(G), then any TDC consisting of r closed trails provides a partition of q 
into its trails and hence (q, P) contributes to (uj r ). On the other hand, if (q, P) contributes 
to (uj r ), then letting G be the graph with the vertices from {1, . . . , N} and the edges given 
by P, we get that G is nimple since q consists of edge-disjoint closed trails. Furthermore, 
G has a TDC consisting of r closed trails and hence (q, P) G c(G). □ 

Proposition 3.4. If c(G) fl c(G') ^ 0, then G is isomorphic to G' . Moreover, if G is 
isomorphic to G' , then c(G) = c(G'). 

Proof. If (g, P) G c{G)C\c{G') , then the construction of q induces a function between the sets 
of vertices of G and G', and P gives the edges of both G, G'. Hence they are isomorphic. 
The second part is true since the definition of c(G) does not depend on 'names' of the 
vertices. □ 

As a consequence we obtain the following. 

Theorem 6. Let uo T be the function defined in / TiOj) . Then 



where the sum is over all isomorphism classes [G] of nimple graphs with at most N vertices 
that have a TDC consisting of r closed trails. 

Next we need to extend the notion of the embedding. A pseudosurface is obtained 
from a surface by identifying finitely many (not necessarily disjoint) pairs of vertices. The 
identified vertices are called singularities of the pseudosurface. A pseudoembedding of a 
graph G is defined to be an embedding on a pseudosurface S such that each singularity of 
S has a vertex of G embedded in it. 
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Proposition 3.5. Let G be a nimple graph. There is a TDC of G consisting of r closed 
trails iff there is a pseudoembedding of G with r faces. 

Proof. The first implication: we construct, from each trail w of TDC, a planar disc bounded 
by a polygon with \w\ edges. Let the boundary be directed according to the direction of 
w. Then we glue the discs together by identifying the oppositely oriented edges with the 
same endvertices. The result is a pseudosurface. 

The other implication: given a pseudoembedding, orient the faces so that each edge 
appears twice, and with the opposite orientation. The orientations of the faces clearly 
form a TDC. 



where the sum is over all isomorphism classes [G] of nimple graphs with at most N vertices 
that have pseudoembeddings with r faces. 

If W is a fat graph, then its dual is the abstract graph W* = (V*,E*) with vertex set 
V* and edge set E*, where V* is the set of (the duals of) the faces of W and E* is the set 
of (the duals of) the edges e of W such that in W*, e = e(fg) has endvertices f,g G V* , 
and in W, e lies on the boundary of faces /, g of W . Note that there is a natural bijection 
between the set of the edges of W* and that of W. 

Next we consider the geometric duals of the pseudoembeddings. 

Definition 3.7. Given a fat graph W and a partition Q of the set of its faces, we denote 
by Gq(W*) the (abstract) graph obtained from W* by the contraction of the classes of Q 
into single vertices. 

A pair (W, Q), where W is a fat graph and Q a partition of the set of its faces, is called 
r-relevant ifW has r vertices, Q has at most N classes, and Gq(W*) is a nimple graph. 

We say that (W,Q) is equivalent to (W',Q'), denoted by (W,Q) ~ (W',Q'), if G Q {W*) 
is isomorphic to Gq>(W'*). 

We call a partition Q of the face set of W trivial and denoted it by if each face of W 
forms a single partition class of Q. 

The following is a consequence of Corollary 13.61 
Corollary 3.8. Let u r be the function defined in pOj) . Then 



□ 



As a consequence of Theorem and Proposition 13.51 we have 



Corollary 3.6. Let uj r be the function defined in UD) . Then 





(11) 



where the sum is over all the equivalence classes [(W, Q)] of the r-relevant pairs where Q 
has a(Q) < N parts. 
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4. Counting planar graphs (Proof of Theorem ED 

In this section we prove Theorem [3] To this end we first derive a relation (Theorem [7j) 
between the enumeration of planar graphs and the Gaussian integration of the generating 
function for edge-disjoint closed trails. We then (in Theorem |H] and Corollary 14. 1 j) derive 
an identity between this Gaussian integration and the Gaussian integration of the ice-type 
partition function for eulerian sets defined in (Tj[|). Theorem [3] follows from Corollary 14.11 

We consider the generating function for edge-disjoint closed trails definded by 

C(M,x,y):=J2 UJ r(M)x r (12) 

r>0 

where u r (M) = E g eW(r) V Ht2 iLeq M e (as defined in (HDD) and W(r) is the set of all 
subsets of edges of D(M) that can be decomposed into r edge-disjoint closed trails. 

Theorem 7. Let s(N) < f be any function in N which approaches to oo slower than N 
when N — > oo. The Gaussian matrix integral of the function ( defined in filfy satisfies the 
following. 




ex HS§^i £l(iV ' r ' e)+ £^ 1) ' £ 

e- N2 (aM,Nzy-\y)) 

' % ' ex p ( E E zT y e ~ r I ^ r ' e ) + E 



s(N) ' 




|Aut(G r 

r>0 e>0 \ [G r ,e] 

where the sum is over all isomorphism classes [G>,e] of connected simple graphs which have 
planar embeddings with r faces, e edges and < N vertices, and E\(N, r, e), E2(N, r, e) are 
of type 0(N~ r ). 

The following theorem formulates the enumeration of planar graphs as the Gaussian 
integration of an ice-type partition function via eulerian sets. 

Theorem 8. Let ( be the function defined in ( TJ1| ) and consider a function r\ which maps 

M eH N to 

r)(M,x,y):= ]T U{A ) x)y^ ,2 \\_M e) (13) 

AcNxN eulerian e£A 

where U{A,x) = He component o/^EE'/M- Then 

(0 = fa>. 

Proof. This is true from the Euler theorem without (•) when we restrict the eulerian sub- 
sets of D to the subsets consisting of pairs of oppositely directed edges with the same 
endvertices. In the integration the remaining terms of both £, t] disappear. □ 
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Note that the function r](M, x, y) defined in (TT3]) is equal to 

r\C\/ 

x — 1 



AcNxN eulerian egA C: component of A 

Thus we have 

(Nz)M/ 2+1 -(Nz)yW 



nMNzy-^y) = Yl (II Me ) 



Nz - v 

AcNxN eulerian eGA C: component of A 

as defined in (pQ). As a consequence of this and Theorem [8] we have the following. 
Corollary 4.1. Theorem\2\is true when we replace ((M, Nzy~ l ,y) by r)(M, Nzy -1 ,y) . 

In the rest of this section we will prove Theorem [7J and then Theorem [3J 
Proof, (of Theorem [7J From Corollary 13.81 we get 

(aM,Nzy-\y)) 
= Y,^)(Nzy- l y 

^ ^ \Aut(G Q (W*))\ 



r>0 [(W,Q)] 



r y F -(w) r W) TT 1-—) 11 (14) 

r>0 [(W,Q)] l<t<a(Q)-l 



NJ \Aut(G Q (W*))\ 



where the sum is over all equivalence classes [{W, Q)} of the r-relevant pairs where Q has 
a(Q) < N parts. 

To see how (1141) gets non-zero contributions from equivalence classes of the r-relevant 
pairs, observe the following. Each equivalence class [(W, Q)} can be represented by a fat 
graph W with some number f(W) of faces satisfying a(Q) < f(W) < 2e(W). 

For further computation we take any representation (W, Q)r of each [{W, Q)] and get 

(((M,Nzy-\y)) 

(■ \ r e(W)—r 

i m z y 

NJ \Aut(G Q (W*))\' 1 j 

We will bound (fl"5l) from above and from below. 
Let Si(N, y, z) and S 2 (N, y, z) be defined as follows: 

2s(iV)\ Q(Q) z'y (H0 ~ r 



SWT, ]T 7V—W+/Wa/«(Q)-/W fl - H^OY 

r e ( W ) — r 



r>0 (w,g)i 



|Aut(G Q (W*))| 

|Aut(G g (W))| 
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where s(N) is any function in iV which approaches to oo slower than N when N — > oo. 
It is straightforward (we recall that that a(Q) < 2e(W)) that 

Si(N,y,z) ' < S{N) ' (C(M,Nzy-\y)) ' < ' S 2 (N,y,z). 

First we study S 2 (N, y, z). Recall that for each representation (W, Q)r, W is a fat graph 
with r labelled vertices and f(W) faces. Define the set S(w,q) r of siblings (W, Q)r of 
each (W, Q)r, where W is obtained from W by relabelling of its vertex labels. Since 
Ew>es (w , Q)R 1 = r!/|Aut(W0| (and notice that v{W) = v(W) = r, e(W) = e(W% f(W) = 
f(W), |Aut(W)| = |Aut(W)| and \Aut(G Q (W*))\ = \Aut(G Q (W'*))\), we have that 

^w'es (WtQ)R \ e(w)+f(w) N a( Q )-f(w) z y 

rl/|A»t(W)l ) |Ant(G,(H~))| 

= V V V N r-e(W')+f(W) N a(Q)-f(W) |Allt(W)| Z r V < W >~ r 

^ ^ ^ \Aut(G Q (W'*))\ r! ' V ; 



Let {Gi, • • • ,G C } be the set of the connected components of Gq(W*). This induces a 
partition of the components of W into W\, . . . , W c , and a partition of the classes of Q into 
Qi, . . . , Q c . Each Qi is a partition of the faces of Wi. 

Notice in dTTJ) that each of the factors, JV r - e ( t O+/0njV a «WC*n = jv«(w')-eflO+a(Q) s 
|Aut(W% |Aut(G Q (W*))|, 2^ e ( w/ ')-'- = ^(W) y e(w')-<W) is multiplicative under disjoint 
union of (Wi,Qi). We have r!/|Aut(W)| ways to relabelling the vertices of W, which is 
equal to 

r \ r\\ x r 2 ! x ■ • • x r c ! 



ri,r 2j -.. ,r c ; |Aut(Wi)| x |Aut(W 2 )| x •■■ x |Aut(W c )|' 

that is, the number of ways of choosing vertex sets for each Wi and labelling Wi. Thus 
we can use the relation between the generating function for not-necessarily connected 
representations and the exponent of the generating function for connected ones to obtain 

S 2 {N,y,z)=Nexp(T 2 (N,y,z)) (18) 

with 

TofiV „ 7l = V V V N r-e(W')+f(W>) N a-f(W) |Aut(W)| Z r y< W >~ r 

W,V,z) 2^ 2^ 2^ \Ant(G Q (W'*))\ r! 

r>0 (W,Q) R W>eS( WtQ) c 



where the second sum is over all representations (W, Q) C R where Gq(W*) is connected. The 
bound — in (fTBl) is imposed, since in S 2 the total degree in y and z (which equals e(W') = 
z{W)) ^ T implies that f(W) < N. Observe that for any nimple graph, f(W) < 2e(W). 
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Since T,wes (WiQ)0 1 = r\/\Aut(W')\, we have 
T 2 (N,y,z) 

V ^ |Aut(G Q (W*))| ^ r! 

t e(W) — T 

V ^ |Aut(G Q (W*))| 

"2^2^ 2^ |Aut(G Q (W*))|' [ ' 

r>0 e>0 1 V VV 7/1 

where has r vertices, e edges, and < N faces. 
We recall that for each r > 0, 

2 - 2<?(W0 + a(Q) - f(W) = v(W) - e(W) + a(Q) = r - e(G Q (W*)) + v(G Q (W*)) 

is independent of the choice of the representation (W, Q)r. 

If (W, Q) C R has g(W) > 0, or g(W) = and Q is non-trivial, then 

N -2g(W) N a(Q)-f(W) _ 

The contribution of such representations to T 2 (N,y,z) is 

Ar2 EE^ 6 " r ^(A,r,e), (20) 

r>0 e>0 

where 

=e;;, 0) ,,o(a'- 1 )^^, (2d 

where the sum is over all the equivalence classes [(W, Q) c ] of the r-relevant pairs where 
Gq(W*) is connected, W has e edges, and g(W) > or Q non-trivial. 

On the other hand, for the remaining representations (W, Q) C R (that is, g(W) = and 
Q is trivial), 

N -2g(W) N a(Q)-f(W) _ ^ 

Hence the contribution of such representations to T 2 (N,y,z) is 

JV 2 yV/f r V 1 —— . (22) 

L^i L^t y ZZ AutfW 7 "*) 

r>0 e>0 [(W,0) c ] pZanar 1 V n 

We can also sum over the duals in (1221) . The dual of a connected fat graph is connected 
and the dual of a relevant fat graph is nimple. Moreover, there is exactly one nimple graph 
on at most one vertex (with no edge) and at most one simple connected graph with one 
vertex and no edge. Hence (1221) is equal to 

N 2 YY z r y e ~ r Y - l - -r , (23) 
Z^Z^ y Aut (G re ) K ' 

r>0 e>0 [G r ,e] ' V ' 
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where the sum is over all isomorphism classes [G r e ] of connected simple graphs which have 
planar embedding (we recall that Q is trivial) with r faces, e edges and < N vertices. 
Putting (PSD, ([2"0]), and (1251) together, we obtain that 



S 2 (N,y,z) = f exp U 2 £E^T r \E 2 (N,r,e) + £ 1 ) I . (24) 

\ r>0 e>0 \ [G r ,e] ' 1 ^ 1 / / 

The analysis for S\(N, y, z) is analogous as to S 2 (N, ?/, z); we only note that the factor 
^1 — ^p- j is multiplicative. We obtain that 

S 1 (N,y,z) = f exp [N^Y,^ \MN,r,e) + £ — -L_ 

V r>0 e>0 V [G r , e ] 1 1 r ' e 



where 



Eerr , / 

[(W,Q) C ] V 7 y 



2s(iV)V (Q) 1 



N J \Ant(G Q (W*))\ 

2s(AQ\ a(Q) \ 1_ 

J^J j |Aut(G r 



where the sum in the first term is as in E 2 (N, r, e), and the sum in the second term is as 
in (F23D. □ 



Ss(N) 



Proof, (of Theorem [3]) From Corollary 14.11 we have that 
e- 7v2 (r / (M,iV^- 1 ,y)) 

'^' exp (5S^H r ' e)+ £,ra))- (26) 

We first study the leading term J2[G T e] 1 / 1 Aut (G^^e) | in ([26]) . We note that to each 
isomorphism class [G>, e ] of connected simple graphs which have planar embeddings with 
r faces, e edges, and v(G r>e ) < N vertices, there are v(G r>e )\ / |Aut(G>,e)| many labelled 
corresponding ones. From this, together with the Euler relation, v(G r>e ) — e + r = 2, we 
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have that 

EE 2 ^ T E lAutfG )\ 

r>0 e>0 [G r ,e] r ' 

~EE E z r y <G T , e )-2 ^ | Aut ( Gr 



|Aut(G re )| r>(G re )! 

r>0 e>0 G r , e labelled 1 v , yl v r ' ' 

EE E w 

r>0 e>0 G re labelled V r ' e/ ' 



Being sorted by a possible number n of vertices, (1271) is equal to 

# z r y n ~ 2 



ni 

n =l r >o 



where p(n, r) denotes the number of labelled connected simple graphs on n vertices that 
have planar embeddings with r faces. 

Concerning the subleading term E 2 (N,r,e) in (|2T)|) . due to the Euler relation, e — r = 
f(W) + 2g(W) - 2, we have 



^^zV" r ^ 2 (iV,r,e) 



r>0 e>0 



^ E E * r y e ~ r E err 



r>0 e>0 



'[(W,Q)=] (Aut^^* 



WQ)1 |Aut(G (W))|- 
If (125]) is sorted by a possible number n where n = /(W) + 2gf(W), it is equal to 



z r y n-2 
r ) j— : 

n>l r>0 



nl 



^^e 2 (iV,n, 

with 

Eerr . 
[<wm° {N ~ '\Ar*(G Q (W*))\> 

where the sum is over all equivalent classes [(W, Q) c ] of the r-relevant pairs where Gq(W*) 
is connected, W satisfies n = f(W) + 2g(W), and g(W) > or Q non-trivial. Note that 
the sum is finite, since e(W) = n + r — 2. 

The analysis of Ei(N, r, e) and ei(N, n, r) is analogous. □ 
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5. Counting directed cycle double covers (Proof of Theorem [3J 

In this section we prove Theorem HI To this end we first formulate the enumeration of 
graphs with a specified directed cycle double cover (DCDC) as the Gaussian matrix integral 
of the generating function £ for even sets with its specified cycle decomposition (defined 
in ©)• We further show that £ is indeed an Ihara-Selberg type function, via a lemma on 
coin arrangements (Lemma 15.51) . 

First we recall from the introduction that a subset A of directed edges of D is even if A 
can be written as a union of edge-disjoint directed cycles of length bigger than two. We 
further denote by A the set of all pairs (q, K) where q is an even set of edges of D and K is 
a decomposition of q into directed cycles of length at least three. The following observation 
is straightforward. 

Observation 5.1. Let (q,K) £ A and let P be a proper pairing of q. Then the cycles of 
K form a DCDC of a nimple graph with vertices {1, . . . , N} and the edges given by the 
pairing P. 

Definition 5.1. Let G be a finite nimple graph with at most N vertices and let C be a 
DCDC of G. Then let c(G, C) be the set of all triples (q, P, C) so that there is a colouring 
d of the vertices of G by colours {1, . . . , N}, where each vertex gets a different colour, 

q = {(d(x), d(y)), (d(y), d{x))\ {x, y} £ E(G)}, 

and P consists of all the pairs ([(d(x), d{y)), (d(y), d(x))]; {x, y} £ E{G)). 

We remark that each such pair (q, C) in Definition 15. II belongs to A and that 

N(N-1)...(N-\V(G)\ + 1) 
|C(G ' 6)I " |Aut(G,C)| ■ 

Proposition 5.2. For the function £ defined in (Tl|) which maps M £ to 

£(M):= H M » 

(q,C)eA eeg 

a term ((q,C),P) with P a proper pairing of q contributes to (£) if and only if there is a 
nimple graph G with at most N vertices where C is its DCDC and such that (q, P, C) £ 
c(G,C). 

Proof. If (q,P,C) £ c(G,C), then C provides a partition of q into its directed cycles and 
hence ((g, C),P) contributes to (£). On the other hand if ((q,C),P) contributes to (£), 
then letting G be the graph with the vertices from {1, . . . , iV} and the edges given by P we 
get that G is nimple. Moreover q comes with a decomposition C into the directed cycles 
and thus (q, P, C) £ c(G, C). □ 

Proposition 5.3. Ifc(G, C) n c(G", C') ^ 0, then (G,C) is isomorphic to (G\C). More- 
over, if (G,C) is isomorphic to (G',C), then c(G,C) = c(G',C'). 
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Proof. If (g, P, C) G c(G,C) fl c(G',C), then the construction of g induces a function 
between the sets of vertices of G and G', and P gives the edges of both G, G'. Hence there 
is an isomorphism from G to G' which preserves the fixed DCDCs. The second part is true 



where the sum is over all isomorphism classes of pairs (G, C) where G is a nimple graph 
with at most N vertices and C a specified DCDC of G. 

5.1. Calculations. The integral (£) counts all the directed cycle double covers of graphs 
on at most N vertices and hence its calculation is an attractive task which need not be 
hopeless. We show next a curious formula for £(M) which identifies it with an Ihara- 
Selberg-type function (see Theorem [TP]) . 

Construction of digraph D'. We first construct a directed graph D' with the weights 
on the transitions between the edges. 

First we split each vertex of D, i.e., we replicate each vertex v, and then we connect 
two w's by a new edge e(v) and we let all the edges of D entering v enter the initial vertex 
of e(v), and all the edges of D leaving v leave the terminal vertex of e(v). If an edge g\ 
enters v in D then we define the weight of the transition w(gi, e(v)) = M gi . We let all the 
remaining transition be equal to one (see Figured], the first two parts). 

Finally, for each pair gi,g 2 of oppositely directed edges of D, say g\ = (uv),g 2 = (vu) 
we introduce a new vertex v g (with g = {u,v}) and we let both gi,g 2 pass through it; 
equivalently, we subdivide both g 1; g 2 by one vertex and identify this pair of vertices into a 
unique vertex called v g (and thus we have new edges (uv g ), (v g v) from g\ = (uv), and new 
edges (vv g ), (v g u) from g 2 = (vu)) (see Figure [U the last two parts). 




□ 




fj2 



U 



V 




I! 



V 



e(v) 



Figure 1. Construction of e(v) and v g 



We let the weights of the transitions at the vertex v g between g\ and g 2 (i.e., between 
(uv g ) and (v g u) and between (vv g ) and (v g v)) be equal to zero, the transitions along g\ and 
g 2 (i.e., between (uv g ) and (v g v) and between (vv g ) and (v g u)) be equal to one, and the 
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transitions between (v g v) and e(v) be equal to M gi and between (v g u) and e(u) be equal 
to M g2 . See an example in Figure [2J 




Figure 2. An example of the construction of D' from D 

In what follows, the directed closed walk is considered not pointed. We let the weight of 
the directed closed walk be the product of the weights of its transitions. 

Observation 5.2. There is a weight preserving bisection between the set of the directed 
cycles of D of length at least three and a non-zero weight, and the set of the aperiodic 
directed closed walks of D' of a non-zero weight which go through each directed edge and 
through each vertex v g at most once. 

Proof. This follows directly from the construction of D' . □ 

Definition 5.4. We define the rotation number p(w) for each closed walk w of D' with 
a non-zero weight by induction as follows: first order the directed edges of D' , say as 
ax, ... , a m , so that the edges e(v),v G V(D) form the terminal segment. Then 

1. If w is a directed cycle, then we let p(w) — — 1. 

2. Let w go at least twice through a directed edge. Let a be the first such edge in the 
fixed ordering. Hence w is a concatenation of two shorter closed walks wi,W2, both 
containing a. If a ^ e{y) for some v then we let p{w) = p{wi)p{w2) ■ If a = e(v), 
then we let p{w) = 0. 

3. If none of l.,2. applies, w must go through the vertex v g (introduced in the definition 
of D' ) at least twice. Then we again let p(w) = 0. 

Theorem 10. Let £ be the function defined in (TJJ). Then 

aM)=n(i-p(p)w(p)), 

v 

where the product is over all aperiodic directed closed walks p in D' and w(p) denotes the 
weight of p. 

To prove Theorem [10] we will need a curious lemma on coin arrangements stated below. 
It has been introduced by Sherman [29] in the study of 2-dimensional Ising problem. 
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Lemma 5.5 (A lemma on coin arrangements.). Suppose we have a fixed collection of N 
objects of which mi are of one kind, m 2 are of second kind, ■ ■ ■ , and m n are of n-th kind. 
Let bk be the number of exhaustive unordered arrangements of these symbols into k disjoint, 
nonempty, circularly ordered sets such that no two circular orders are the same and none 
are periodic. For example let us have 10 coins of which 3 are pennies, 4 are nickles and 3 
are quarters. Then {{p, n), (n,p), (p, n, n, q, q, q)} is not a correct arrangement since (p, n) 
and (n,p) represent the same circular order. If N > 1 then Y2iLi(~ = 0- 



Proof of Lemma \5.5\ The lemma follows immediately if we expand the LHS of the following 
Witt Identity and collect terms where the sums of the exponents of the z^s are the same. 

Witt Identity (see [T7]): Let Zi,...,z k be commuting variables. Then 
H (1 - = 1 - z x - z 2 - ... - z k , 

mi ,...,mj.>0 

where M(mi, ....,mk) is the number of different nonperiodic sequences of z^s taken with 
respect to circular order. □ 



Proof of Theorem QJJ We first show that the coefficients corresponding to the products of 
variables where at least one M e , e ^ e(v), appears with the exponent greater than one, are 
all equal to zero. 

Let us denote W(p) = —p(p)w(p). Let A\ be the set of all aperiodic closed walks p such 
that di appears in p. Each p G Ai has a unique factorization into words (Wi, Wk) each 
of which starts with a\ and has no other appearance of a\. 

Let S be a monomial summand in the expansion of rLeAi(l + W(p)). Hence S is a 
product of finitely many W(p),p G A\. 

Each p G Ai has a unique factorization into words defined above. Each word may appear 
several times in the factorization of p and also in the factorization of different aperiodic 
directed closed walks. Let B(D') be the set-system of all the words (with repetition) 
appearing in the factorizations of the aperiodic directed closed walks of D' . 

It directly follows from Lemma 15.51 the lemma on coin arrangements, that the sum of 
all monomial summands S in the expansion of n p eAi(l + W(p))> which have the same set 
of 'coins' of B(D') with more than one element, is zero. Hence the monomial summands 
S which survive in the expansion of IIpeAi^ W(p) all have the set of coins of B(D') 
consisting of exactly one word. Such S cannot have a\ with exponent bigger than one. 
Now we can repeat the same consideration for the other edges different from e(v),v G V. 

Hence the only terms of the expansion of the infinite product that survive have all 
M e , e e(v), with the exponent at most one. 

We know from Observation 15.21 that the collections of the edge-disjoint directed cycles 
of length at least three in D correspond to the collections of the aperiodic directed closed 
walks of D' where each edge of D' and each vertex v g appear at most once; by above, these 
exactly have chance to survive in the infinite product. 
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Each term of £(M) may be expressed several times as a product of aperiodic directed 
closed walks of D', but only one such expression survives in the infinite product since if a 
closed walk goes through an edge e(v) or through a vertex v g more than once, its rotation 
is defined to be zero. Hence is counted correctly in the infinite product. □ 

Remark 5.6. Let us write p(p) = (— l) rot ( p ). Without the zero values of p(p), the func- 
tion rot(p) is additive when we 'smoothen' p into directed cycles. The integer lattice 
generated by the directed cycles has a basis which may be constructed e.g., from the 
ear-decomposition [TT]; the function rot(p) may be split into contributions of the edge- 
transitions for the basis, and since it is a basis, it may be split into contributions of the 
edge-transitions also for all the directed cycles. Hence if the additivity property holds, 
rot(p) may be split into the contributions rot(t) of the edge-transitions t for the aperiodic 
closed walks. Hence we would have 

n(i - P (pmp)) = na - n (-ir ( M*))- 

p p t transition of p 

This formula transforms the infinite product into the Ihara-Selberg function. It was studied 
by Bass in pQ who proved that it is equal to a determinant. A combinatorial proof was 
given by Foata and Zeilberger in [TP] . 

Having the zero values of p(p) it is not clear how to split the rotation to individual 
edge-transitions. A determinant-type formula may however exist. 
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